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We propose a learning technique for MIMO secondary users (SU) to spatially coexist with Primary 
Users (PU). By learning the null space of the interference channel to the PU, the SU can utilize idle 



> 

\ degrees of freedom that otherwise would be unused by the PU. This learning process does not require 



any handshake or explicit information exchange between the PU and the SU. The only requirement is 
that the PU broadcasts a periodic beacon that is a function of its noise plus interference power, through 
' a low rate control channel. The learning process is based on energy measurements, independent of 



the transmission schemes of both the PU and SU, i.e. independent of their modulation, coding etc.. 
The proposed learning technique also provides a novel spatial division multiple access mechanism for 



equal-priority MIMO users sharing a common channel that highly increases the spectrum utilization 

■ 

compared to time based or frequency multiple access. 



I. INTRODUCTION 

The emergence of Multiple Input Multiple Output (MIMO) communications opens new di- 
rections and possibilities for spatially sharing wireless channels QjJ— 13[| . Consider a scenario of 
two independent MIMO communication systems that share the same flat fading MIMO channel 
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Fig. 1. Blind spatial division multiple access for MIMO users with equal priority. The matrix Hi t j,i 7^ j £ {1,2} are 
unknown to both users. The objective of the two users is to leam the null space of Hij, i ^ j 6 {1, 2},. 



as depicted in Figure [Q Assuming that each user has more antennas at the transmitter then the 
maximum number of antennas that each one has at the receiver, they can share the channel 
without interfering to each other by using orthogonal spatial dimension. This spatial sharing 
is even more appealing in MIMO Cognitive Radio (CR) networks [|4j-l6D since it enables a CR 
MIMO Secondary User (SU) to transmit a significant amount of power simultaneously as the 
PU without interfering with him by utilizes spatial dimensions that are not used by the PU. This 
spatial separation requires, in both CR and MAC, that the interference channel be known. In the 
MAC (see Fig. Q]), it means that H 2 i and H 12 be known to user 1 and 2 respectively, while in 
the CR case it is sufficient that the SU, say user 2, knows H 2 i. This information can be achieved 
by conventional channel estimation techniques that require a high level of cooperation, including 
handshake, transition of a known synchronous training sequence and the use of matched filters 
for each receiver antenna. In the MAC scenario, this process needs to be applied twice were at 
the first time one of the systems transmits a training sequence while the second estimates the 
channel and transmits the estimation back to the other system, then it is repeated where the two 
systems exchange roles. During these processes, each system must stop its data flow unless it 
is capable of full duplex, i.e. transmitting and receiving simultaneously at the same time and on 
the same frequency band. Although complicated, this channel estimation can be carried out in 
MAC since both users are equal priority. In CR on the other hand, this is far more complicated 
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since nobody expect PUs to stop their reception and perform channel estimation for unlicensed 
secondary users. Thus, acquiring and/or operating without knowing the interference matrix to 



the PU is a major issue of active research [M-lnll. in CR. Note that every solution that is good 



for CR problem can be utilized to the MAC problem. Henceforth, we consider the problem of 
interference channel learning in the context of CR. 

We consider the underlay CR paradigm ifh?!. that is, the SUs are constrained not to exceed a 
maximum interference level at the PU. The optimal power allocation for the case of a single SU 
who knows the matrix H 2 i in addition to its own Channel State Information (CSI) was derived 

n A 

by Zhang and Liang Q4Q . In the case of multiple SUs, Scutari at al. [6] formulated a competitive 



game between the secondary users. Assuming that the interference matrix to the PU is known 
by each SU, they derived conditions for the existence and uniqueness of a Nash Equilibrium 
point to the game. Zhang et al. [|80 were the first to take into consideration the fact that the 
interference matrix Hi 2 may not be perfectly known (but is partially known) to the SU. They 
proposed Robust Beamforming to assure compliance with the interference constraint of the PU 
while maximizing the SU's throughput. Another work for the case of an unknown interference 
channel with known probability distribution is due to Zhang and So |90 who optimized the SU 
throughput under a constraint on the probability that the interference to the PU be above a given 
threshold. 

A very appealing solution concept for CR in general and MIMO CR in particular, is that the 



indly without a handshake and without 
Proposed such a solution in the case 



SU would be able to mitigate the interference to the PU b 
using conventional channel estimation techniques. Yi Jil l 
where there is a channel reciprocity between the PU's transmitter and receiver in which the SU 
listens to the PU signal and estimates H 12 's null space from its second order statistics. This 



work was enhanced by Chen et al. in II 1011 . Both works require channel reciprocity and therefore 
are restricted to a PU that uses Time Division Duplexing (TDD) . Once the SU obtains the null 
space of H 12 it can transmit within this null space without interfering with the PU. 

Beside the channel reciprocity case, obtaining the value of by the SUs (i.e. the interference 
channel to the PU) requires the PU to participate in the SU's estimation task. This task requires 
that the SU transmits a training sequence, from which the PU estimates and feeds it back to 
the SU. Such cooperation increases system complexity overhead, since it requires a handshake 
between both systems and in addition, the PU needs to be synchronized to the SU's training 
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Fig. 2. The addressed cognitive radio scheme. The matrix H12 is unknown to the secondary transmitter and vi(fe) is a 
stationary noise (which may include stationary interferences) . The interference from the SU, Hi2X2(fe), is treated by the PU as 
noise, i.e. no interfere cancellation is performed. The SU obtains a closed form expression for the null space of the interference 
channel to the PU H12 by measuring the variation of qik) resulting of finite set of transmitted signals {x2{k)}J =1 . 

sequence. This required cooperation is one of the major obstacles to deployment of MIMO CR 
systems. 

The objective of this paper is to design a simple procedure based on minimal cooperation by 
the PU such that a MIMO SU will be able to meet its interference constraint without explicitly 
estimating the matrix and without burdening the PU with any handshake, estimation and/or 
synchronization associated with SUs. Consider the problem depicted in Fig. |2] In this scheme 
the PU, although active, is not necessarily aware of the SU. Its role in the SU's learning process 
is limited to broadcasting a single one-dimensional beacon through a low rate control channel. 
This beacon is a function of the PU's noise plus interference. The advantage of this technique 
over conventional channel estimation techniques is that it does not require a handshake and 
synchronization between the secondary and the PUs and can be implemented using only energy 
measurements. This is also a very appealing property for interference mitigation between two 
MIMO users (i.e. "multiple access") since it makes the information exchange mechanism between 
the two users that is needed for them to share the same channel very simple. 

The remainder of this paper is organized as follows: Section HI formulates the problem. Section 
lITll presents the Energy Based Cannel Learning (EBCL) algorithm for interference mitigation in 
the primary -secondary user CR scenario. Section |IV]discusses the implementation of the EBSL 
algorithm in spatial channel sharing between two independent MIMO users of equal priority. 
Section |V] presets numerical results. 
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II. Problem Formulation 

Consider a flat fading MIMO interference channel with a single PU and a single SU without 
interference cancellation, i.e. each system treats the other system's signal as noise. User's i's 
i G {1,2} received signal is given by 

yi{k) = KaXi(k) + HijXj-(fc) + v^k), fceN (1) 

where j G {1,2}, j ^ i, H iq G C riXtq and Vj(A;) is a zero mean stationary noise. In this 
paper all vectors are column vectors. Let A be an / x v complex matrix, then, its null space 
is defined as A/"(A) = {y G C v : Ay = 0} where = [0, ...,0] T and its column space 
C(A) = span(a 1; . . . , a^) C C l . We assume that user 1 is the PU. The secondary user (user 2) 
is allowed to transmit as long as the interference does not exceed a maximum level at the PU, 
i.e. 

||H 12 x 2 (A;)|| 2 <77, (2) 

where rj = represents the case where the SUs are allowed to transmit only in the null space 
of the matrix H 12 . 

Since the secondary user is MIMO it can share the channel without interfering with the PU if 
it uses spatially orthogonal degrees of freedom. In particular, the SU will not interfere with the 
PU if its transmitted signal x 2 satisfies x 2 G A/"(H 12 ). The main obstacle in using this technique 
is that it requires knowledge of A/"(Hi 2 ). The matrix Hu is known only to the PU, and the 
matrix H 12 is unknown to both the PU and the SU; hence its estimation requires cooperation 
between the two users. The state of the art in MIMO channel estimation techniques requires 
that the SU transmits a training signal that is known to the PU. The PU then estimates the 
channel using a matched filter. Other techniques that are not based on a known deterministic 



signal waveform are the blind MIMO channel estimation techniques H13M15I. e.g.] in which the 
receiver uses the received signal statistics, i.e covariance matrices and higher order comulant 
tensors, to estimate the channel. These approaches require an extensive set of measurements and 
processing at the receiver side (the PU's receiver in this case). After the PU obtains an estimate 
of Hi 2 he transmits it to the SU. This kind of "service" provided by the PU to the SU is highly 
undesirable due to the overhead and cooperation required on the part of the PU. Thus, reducing 
the role of the PU in this channel learning phase will make CR technology more attractive for 
practical applications. 
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Our objective is to derive a simple procedure for the SU to learn the null space of the matrix 
H 12 such that the PU would not need matched filters or to make extra measurements other than 
those required for its usual operation. We would also like to reduce the amount of processing 
at the PU and above all we would like the SU to obtain the null space of H 12 without having 
a handshake with the PU and even without the PU being aware of the SU. We denote 

G = H* 2 H 12 e C t2Xt2 (3) 

and divide time into iV-lengfh intervals referred to as transmission cycles. In each transmission 
cycle, the SU transmits a constant signal (this is required only during the learning process) , i.e. 

x 2 ((n-l)iV) = x 2 ((n-l)7V+l) 

= • • • = x 2 (iVn - 1) = x 2 (n) 
while the PU measures its total noise plus interference. It then broadcasts to all of the users in 
its vicinity the one dimensional signal q(n) that satisfies the following assumption. 

Assumption 1: There exist some K E N such that the value of ||H 12 x 2 (n)|| 2 can be extracted 
up to an arbitrary scalar factor a > from {g(/)}" =0 , for every 1 < n < K. 

Note that from the SU point of view, knowing H i2 at transmitter is equivalent to knowing G, 
which is defined in ©. The problem of learning the G from {q , (0}™=0' referred to as the energy 
based cannel learning problem, is depicted in Figure |3] Note that as long as a is constant for 
every 1 < n < K, the function q{n) can be measured via energy detectors since a is arbitrary. 

A natural choice for a beacon that satisfies Assumption Q] is the following: 

1 Nn 

= w Yl ElHy^-Hnx^)!! 2 } (5) 

k=(n-l)N+l 

where x^/c) is the decoded signal. This beacon is transmitted at time instances k = nN, n G N. 
If we neglects the decoding errors, (i.e. iti(k) = x-i(k)) we obtain 

Qin) = ^£fj;;V E {l|H 12 x 2 H +v(fc)|| 2 ) 

= ||H 12 x 2 (n)|| 2 + Tr(E{ Vl (£;K(A;)}) (6) 
= x 2 (n)Gx 2 (r2) + c 

We will now show that this beacon satisfies Assumption [H i.e. that the secondary user can 
extract a||Hi 2 x 2 (n)|| 2 from {g(/)}" =0 . This is done as follows: At the beginning of the learning 
process (n = 0) the SU transmits x 2 (0) = 0, that is, it does not transmit. Let a > be the 
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Fig. 3. Block Diagram of the energy based cannel learning problem. The SU objective is to learn the null space of H12 by 
inserting a series of {x2(n)} ng N and measuring the output q(n). The only information that can be extracted by the SU is that 

||Hi 2 x 2 (n)|| 2 > ||Hi 2 x 2 (OI| 2 ifq(n) > q(l) for every (k - 1)K < l,n< kK where k G N. 



magnitude of the control channel from the PU to the SU. Then, at time k = the SU measures 
aq(0) where q(0) = Tr(E{vi(A;)v*(A;)}. For n > 0, the SU transmits the signal ic 2 (n) and at 
time k = nN it measures the aq(n) broadcast by the PU. The SU then obtains a||Hi 2 X2(n)|| 2 
by subtracting aq(0) from aq(n). Note that a may be unknown to the SU and that the only 
requirement is that it be constant during the learning process. 
In practice, the beacon will be based on the sample average 

1 Nn-l 

^ n ) = iV E l|yi(*)-Hu*i(A0|| 2 (7) 

k=(n-l)N 

which depends on the averaged value of ||z(A;)|| at the nth cycle where 

z(fc) = H 12 x 2 (A;) + v 1 (A;) (8) 

It is important to stress that the function q(n) is calculated entirely from yi(fc). Therefore it 
is calculated by the PU processing unit after decoding its signal x^/c) without any additional 
measurements. 

In the next section we will show how the SU can learn the null space of the matrix H i2 from 
the measurements {q(n)} t2 =1 , where t 2 is the SU's number of transmit antennas. 

III. The energy based cannel learning algorithm 

In order to obtain Hi 2 's null space it is sufficient to calculate G's null space (where G is de- 
fined in Q}). The following proposition expresses the matrix G as a function of {x 2 (n)Gx 2 (n)} n 2 = 
where each x 2 (n) is a different transmitted signal. 
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Proposition 1: Let 5(A, x) = x*Ax and m (6 l , 0) be a t 2 — dimensional column vector whose 
entries are all equal to zero except of the /th and mth entry, which are equal to cos(6 l ) and 
sin(6>), respectively, i.e. 

r,, m (0,^) = [O,-..,O,cos(0) J 

0,-- - ,O,e-^sin(0),O, ••• ,0] T 

The entries {gi, m }\ 2 m= i of the matrix G = H^ 2 H 12 are given by 

g u = S(G,n, m (0,0)) (10) 
»W = Q, m (7r/4,0) (11) 
3(»,m) = -Q, m (7r/4,7r/2) (12) 

where 

ci, m (9, 4>) = (gi,i cos 2 (9) + # mjm sin 2 (6>)) 

-SiG^Mfl) (13) 

Proo/- Note that 

^(G, r;, m (0, 0)) = cos 2 (6>) sin 2 (6») |y m>m | 

- |^ >m | sin(20) cos(0 + Z^, m ) 
from which (flOl) follows. By substituting ([14)) into (fl"3l) we obtain 

cj,m(0, 0) = sin(20) |^, m | cos (0 + Z^, m ) (15) 

from which £□) and (JT2J follow. □ 
The EBCL algorithm provides a closed form expression for the matrix G. For every x, the 
value of ||Hx|| 2 can be obtained by transmitting x(n), receiving q(n) and subtracting q(0) from 
it, i.e. 

||HxH|| 2 = gH-g(0) (16) 

From Proposition CD it follows that the matrix G can be obtained precisely by t\ transmission 
cycles. The CF-BNSL algorithm is described in Table U After obtaining the matrix G, its null 
space can be calculated offline at the secondary transmitter's processing unit. Once the SU knows 
the null space of the interference channel to the PU's transmitter it can transmit freely as long 
as its transmitted signal is restricted to its null space, i.e. x 2 G A/"(H 12 ). 
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TABLE I 
The EBCL algorithm 

function G=EBCL 

Set b = S(G,0); 
for / = 2, ...,t 

Set flu = S(G, 1^(0,0)) -6; 

end for 

for I = 1, ...,t- 1 

for m = / + 1, t 

Set a;, m = S(G,r ; , m (7r/4,0) - b; 
Set A, m = S(G, r,, m (7r/4, tt/2)) - 6; 
Set ci = Calc_c(gu,g mm ,ai,n/A); 
Set c 2 = Calc_c(gii,ffmm,/3;,7r/4); 

end for 

end for 
end EBCL 

function: c =Calc_c (gi,g2,a,8) 

c = gi cos 2 (8) + g-z sin 2 (#) — a; 
end Calc c 

The advantage of the proposed scheme (see Fig. |2]) is that the PU, although active, does not 
have to be aware of the SU. Its role in the SU's learning process is limited to broadcasting 
periodically the beacon q(n) through a low rate control channel to all of the secondary users 
in its vicinity. Thus, in order to implement the EBCL algorithm, the secondary user needs only 
to detect and measure q(n)'s energy in every transmission cycle without having a handshake 
with the PU. Recall that the only condition required for the EBCL is that Assumption Q] holds. 
This assumption holds even if there are multiple secondary users in the system as long as 
their interference to the PU is stationary. Thus a new secondary user can join the network while 
multiple SUs coexist with the PU in a steady state, i.e. they are not varying the spatial orientation 
or their transmitted signal. 
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IV. EBCL ALGORITHM FOR SPATIAL DIVISION MULTIPLE ACCESS 



The fact that the CF-BNSL algorithm is based entirely on energy measurement and not on 
matched filters makes it very appealing for implementation as a blind spatial division multiple 
access technique for MIMO users with equal priority (see Figure 3), that is, © is no longer 
required. This simplifies the coordination between the two users as follows: At the first stage, 
there is a handshake between the two systems in which it is decided which system begins with 
learning and which provides feedback. Assume that system 2 begins with learning while system 
No. 1 feeds back its measurements. Then system No. 2 transmits a signal x 2 (n) while system 
No. 1 measures and feeds back the beacon in ©. This way, system 2 learns the matrix Gi by 
applying the CF-BNSL algorithm. This process is then repeated where both systems exchange 
roles such that system 1 learns G2. Thus, if system 1 and 2 restrict their transmission to A/"(H 2 i) 
and A/"(Hi 2 ) respectively, they do not interfere with each other and create in effect a Spatial 
Channel Sharing (SCS) . 

An important question that arises is whether the spatial channel sharing is worth the effort 
of null space learning. Recall that in the primary- secondary user CR scenario the SU must be 
invisible to the PU. This fact makes the learning of jV(H 12 ) worthwhile because, as long as the 
channel remains unchanged, the SU is operating freely without colliding. This is not the case for 
MIMO users of equal priority. They can choose not to mitigate interference at all or to share the 
channel using a much simpler multiple access scheme such as Frequency Division Duplexing 
(FDD), which is static and does not require null space learning. In the sequel it is shown that 
the SCS provides a much better spectrum utilization (in terms of degrees of freedom) than FDD 
if both systems have a sufficient number of antennas at the transmitter. 

In the sequel it is assumed that 1 < % ^ j < 2, tj > r, and that the EBCL algorithm is 
performed by both users. Let 



(17) 



and let 



W,A t W* = G t 



(18) 



be its eigenvalue decomposition. Then user i's pre-coding matrix Tj is given by 



Tt = [w, 



(19) 
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where is Wj's gth column and q x , q 2 , qti- rj are the indexes that chose the eigenvectors 
that correspond to G/s Null space, i.e. 

w;G,w 9l = • • ■ = w*_ r .GjW tj ^ = (20) 

The following proposition shows that for Zero-Mean Spatially Whitej (ZMSW) channels that 
satisfy tj > r J? the EBCL results in a free interference r, x (tj — r^-ZMSW channel for each 
user. 

Proposition 2: Assume that H iq , q,i E {1, 2} are rj x t q (ZMSW) channels that are indepen- 
dent of each other and satisfy tj > rj. Let tin be user z's equivalent channel when both users 
apply the CF-BNDL algorithm i.e. = HjjTj where T, is users z's pre-coding matrix defined 
in dT9l). Then, H& is an r, x (U - rf) ZMSW channel. 
Proof: See Appendix lAl 

Proposition |2] implies that if t 4 > rj + rj for i ^ j G {1, 2}, the difference between SCS using 
the EBCL algorithm compared to the case where there is no interference is equivalent to not using 
rj antennas. Furthermore, both users would not lose degrees of freedom compared to the case 
where there is no interference since rank(Hjj) = r, a.s.,and rank(Hjj) = min{r,;,ij — rj} a.s. 
which are equal if U — r, > r 4 . The following theorem extends the last statement for a wider 
range of channel types. 

Theorem 3: Assume that H iq i,q E {1,2} are independent (i.e. independent of each other) 
random matrices defined on the same probability space (Q,J-, P) such that vec(Hjj) is a 
continues random vectoi^l for i = 1,2. Let cL = rank(Hjj) be user i's number of degrees 
of freedoms if he is operating alone, and let d™ be user i's number of degrees of freedom when 
both users apply the EBCL algorithm, i.e. d" = rank(Hjj), where = HjjTj and T, is users 
i 's pre-coding matrix defined in (fT9l ). Then, dj = df a.s. if U > + rj. 
Proof: See Appendix [Bj 

'it means the the entries of the matrix H are i.i.d. zero-mean unit- variance circular Gaussian random variables [see e.g. [TfiL 
Section 10.1]. 

2 A ^-dimensional complex random vector x is said to be continuous if it can be written as x = xr c + ixi m where x = 
[x^XjnJ such that x is a continuous 2t-dimensional random vector, i.e. x has a probability density function with respect to 
the Lebesgue measure. 
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Fig. 4. Graphical illustration of the space Ai. Assuming that all matrices are full rank and that the secondary user has more 
antennas ate the receiver than at the transmitter, i.e. n > ti, then C x (Hn) 7^ 0. Then, the subspace A2 that maps signals to 
C x (Hn) can be used by the SU without interfering with the PU. A necessary and sufficient condition is that t-z > ti. 



V. Obtaining Additional Degrees of Freedom 

Constraining the SU to transmit only in A/"(H 12 ) may be inefficient in some cases. Consider 
a scenario where the PU has more antennas at the receiver than at the transmitter i.e. ti < 
r\ = rank(Hn) and full CSIR of its own channel Hn. Then, the PU's signal of interest at 
the receiver, that is HnXi, can lie only in the ri-dimensional subspace C(Hn) C C n . This 
redundancy can be utilized by the SU to obtain additional degrees of freedom by transmitting 
x 2 G -/V(Hi 2 ) + aE where A 2 = {x 2 G C* 2 : H 12 x G ^(Hn)} (see Fig. g] for illustration). 
If all matrices are full rank, a necessary and sufficient condition for A/"(Hi 2 ) + A2 7^ is that 
r 2 > n. Note that the subspace J\f(B. 12 )+A 2 is equal to A/"(P Hll Hi 2 ) where P B = B*(BB*) f B 
is the projection matrix into the column space of B (which is equal the rang of B) and (-)t 
represents the pseudo inverse operation. These extra degrees of freedom can be obtained by the 
EBCL algorithm with no additional cost. The only modification required is for the PU to project 
z(t), defined in ([8]), into C(H n ) while the rest of the algorithm remains the same, i.e. to replace 
z(t) = yi(t) — xi(t) with z(t) = Pn 11 z(t). This idea can also be implemented in the case of 
two users with equal priority that is described in Section ITVl 

3 The sum of two vector subspaces is the vector space created by the sum of all the vectors in these two subspaces, i.e. let B 
be a vector space and let Si, B2 be two vector subspaces of B, then Bi + B2 = {x G B : x = y + z, y £ Bi, z = £ B2}. 
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VI. Numerical Examples 

To determine the value of null space learning in this setting we turn to simulations. Figure |5] 
compares the rate gain of SCS over that of FDD in a two-user symmetric MIMO interference 
channel without interference cancellation. By symmetric we mean that ti = t 2 , T\ = r 2 and 



that Hjj,z = 1,2 are ZMSW channels as well as H^,z ^ j, G {1,2}. Figure [5(a)1 shows that 
for t = 4 and r = 2 , the SCS outperforms the FDD, i.e. the SCS's rate gain is higher than 
that of the FDD. Furthermore, in the high SNR regime the SCS rate converges to the channel 
capacity without interference, i.e. the rate of a single user occupying the entire channel, as long 
as t > 2r, as shown in Figure |5(b)| From this we conclude that in the FDD scheme, each user 
exploits only half of its degrees of freedom, whereas in the SCS scheme both users exploit all 
of their degrees of freedom (as long as > r, +rj) and the only performance loss is due to the 
restriction of the transmit signal to J\f(Hji). 

It is important to stress that knowing G can be utilized for a more sophisticated channel 
sharing than the SCS. For example, suppose that in addition to transmitting in ^(Ph^H^), 
system i wishes to use also part of its orthogonal compliment A/'" L (PH 33 H?j). This of course 
creates interference to system j. However by choosing eigenvectors that correspond to G's 
lowest eigenvalues, system i can balance between its performance gain and the interference to 
system j. To show that explicitly, let VSV* be the eigenvalue decomposition of G, where 
S is a real nonnegative diagonal matrix that contains G eigenvalues in decreasing order, i.e. 
ci > °2 > • • ■ > cfd > where d < U. Then the eigenvector that corresponds to a a (i.e. 
V d 's column) produces minimum interference to system j. This way, system i can balance 
between choosing eigenvectors that provide it with the best performance gain and minimizing 
the interference to system j. 

VII. CONCLUSIONS 

We proposed a blind technique for MIMO SUs to spatially coexist with PUs based on minimal 
cooperation from the PU. This cooperation does not require additional sensing by the PU and 
is carried out by calculating the power of the PU's total noise plus interference. This value is 
broadcast via a low rate control channel to all of the SUs in its vicinity (beacon). By doing so, 
the PU enables the SU to utilize unused degrees of freedom. 

The advantages of the proposed technique are: 
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Fig. 5. Comparison between blind spatial division and FDD/TDD in a symmetric MIMO interference interaction. The matrices 
Hi,q i,Q € {1, 2} are i.i.d. complex Gaussian with zero mean and unit variance. Both user's number of received antennas is 2. 
The interference expected power is 10.5dB lower than the expected signal power for both users. The vertical axis represents the 
ratio between the achievable rate to the rate obtained via uniform power allocation over the entire band/time. In Subfigure (a) 
The horizontal axis represents expected received SNR while the number of transmit antennas for each user is 2. In Subfigure 
(b) the horizontal axis represents the number of antennas at the transmitter while the expected power at each receiver is 140 dB. 



1) The SU operates autonomously and independently of the PU (as long as the PU transmits 
the defined beacon). 

2) The PU produces the beacon from information that already exists in all communication 
systems, i.e. from the PU's decoded signal and its received signal. 

3) The entire learning process is based on energy measurements, independent of the transmis- 
sion schemes of both the PU and SU, i.e. independent of their modulation, coding etc.. 
This flexibility is very important in CR networks which are inherently ad-hoc. 

4) The entire learning process takes t 2 2 transmission cycles where t 2 is the number of the 
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SU's transmit antennas. 

5) The proposed technique is easily applicable to CR networks with one PU and multiple SUs 
as long as only one SU performs the learning procedure at a time while the other SUs don't 
change their spatial power allocation. In practice, this is not a problem since the learning 
process takes only t\ transmission cycles. 
For the same reasons the proposed scheme can be easily implemented for spatial channel 
sharing of two independent MIMO secondary users of equal priority. We demonstrated that if 
both users share the channel using the CF-BNSL algorithm: 

1) They don't loss degrees of freedom while gaining an interference free MIMO channel. 

2) In case of for zero-mean spatially-withe Gaussian channels and U > rj, then the SCS results 
in a free interference r; t x (ti — r,) -zero-mean spatially-withe Gaussian channel for each 
user. 



Appendix A 
Proof of Proposition [2] 

Without loss of generality we set i = 1, j = 2 and denote Hi = HnWi. Since H 2 i is ZMSW 
channel, the random matrix G 2 (defined in (fT7l)). by definition, is a central Wishart Matrix. Thus, 



Wi (defined in (1181) ) is a unitary matrix that is uniformly distributed over the manifold of unitary 
matrices in C* lXtl [see e.g., [l7|, Lemma 2.6]. Since the channel H n is ZMSW it is bi-unitary 
invariant Ill7ll . that is UH n V's distribution is unchanged for any unitary matrices U, V. Thus, 
for every Wi, the conditional distribution of Hi's is equal to Hn, i.e. P(Hi/Wi) = P(Hu). 
Therefore, given W x , Hn entries are i.i.d. zero-mean unit-variance complex Gaussian random 
variables (i.e. ZMSW channel) and because this distribution is not a function of W 1; the marginal 
distribution of Hi is the same, i.e. P(Hi) = P(Hi/Wi). It follows that Hi is a n x t\ ZMSW 
channel and therefore Hi (which is composed of some t\ —r 2 columns of Hi) is an r x x (i x — r 2 ) 
ZMSW channel. 



Appendix B 
Proof of Proposition [3] 

In this proof we shall use some special notation. Matrices will be denoted by italic upper 
case letters (i.e. the channels H iq ,i,q = 1,2 are now denoted by H ig ,i,q = 1,2) while 
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random matrices will be denoted by boldface upper case letters. We will make not notational 
distinction between scalars and vectors and denote both with lower case italic letters. Random 
vectors/variables will be denoted by boldface lowercase letters. Without loss of generality, we set 
2 = 1 and denote Hj^ = H and H^ 2 = H. Let h q , h q be PT s and H's, qth columns respectively 
and H_ g be the t\ x (r x — 1) matrix that results from deleting H 's gth column. 

The Theorem is first proven for real matrices. In this case h q , h q : f2 — > R' 1 are Borel 
measurable functions. If r\ < t\ — r 2 , user 1 losses at least one degree of freedom iff there exists 
a sequence of scalars {a g } ? Li not all zero such that Y^q=i a q^ q e N ± (H 2 i) = span(/i 1 , h r2 ). 
The later is equivalent to the following statement: There exists 1 < q < r\ such that h q G C(B_ g ) 
where B_ f/ = [H_ g ,H]. Using the sub-additivity of measures 

P(df < dO < P ( Q h ? G C(B_ g ) ) < P (K e C(B_,)) (21) 



13 



Note thal_ 

P(h g G C(B,)) = / P(h, G C(B_ g )|H)dP(a;) (22) 

It remains to show that P(h q G C(B 5 )|H) = 0, a.s. By hypothesis, H is independent of H, 
thus P(h q G C(B 9 )|H) = P(h q G C(B q )),a.s. Now recall that P H is absolutely continuous 
with respect to the Lebesgue measure, that is, Ph << m' iri where m k is the k- dimensional 
Lebegsue measure. Let Q(Z) = {[x,Y] : x G C([Y,Z]),Y G R*^^" 1 ), Z G R* lXr ' 2 } and let 
Q Y (Z) = {x : [x,Y] G Q(Z)} be Q(Z)'s F-section. Then for every Z G M* lXr2 

m ilXri (Q(Z)) = / m tl (Q Y (Z))dm tl ^ ri ~ 1 \Y) (23) 

[see e.g. LL9|, Theorem 2.36] and since for every Z, Y, Qy(Z) is a vector subspace of R* 1 whose 
dimension is at most r\ + r 2 — 1 it satisfies m' 1 (Q Y (Z)) = (recall that r 2 + T\ < ti). This 
establishes the desired result for real channel matrices. 

To extend this result to complex matrices, note that h q = h_ Re +ih ^ and h q = h q Re + ih q lm 
where h q Re , h Im , h Re , h Im : H — > R* 1 are Borel measurable functions. Furthermore, the 
vector space C* 1 is isomorphic to R 2 ' 1 , that is, there exists a bijective mapping (one to one and on 

4 The existence of a conditional probability measure P(-\b. g )(ui) for each lj 6 is due to the fact that all random vectors are 
assumed to be E' 1 -Borel measurable. Such probability measure is termed regular conditional probability [see e.g. fisll . 
5 Let fi, v be two measure defined on the same measurable space (X, M) , then /i << v if v(A) = =>■ fJ-(A) — 0. 



Febraary 3, 2013 



DRAFT 



17 



to) from one to the other which in this case is given by ip(x) = [Re (or), Im(x )] where x G C* 1 . 
Let = [ — Im(x T ), Rm(i T )] T then C(B q ) is mapped into V q = span('0(/i 1 ), V>(7ii), 

tl}(h r2 ),ip{h r2 ), ipih), ipih), V'(^-i), ip{h q -i), ip(h q+ i), tp(h q+1 ), V(^n), ip{hri))- Thus > ^ g G 
C(£ g ) iff V(h 9 ) G V or ^(h* 1 ) G VQ Because ^{h} 1 ) and ^(h- 1 ) are orthogonal, h q G C(B q ) is 
equivalent to ip(h ) G V 1 - or ^(h ) G V. Henceforth the proof is identical to the real case since 
m 2tl (y) = m^iy-*-) = and because H is a continuous random matrix. 
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